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ABSTRACT 



^ , In order to probe the activity of the inner disk flow and its effect on the neu- 

Q> I tron star surface emissions, we carried out the timing analysis of the Rossi X-Ray 



Timing Explorer observations of the millisecond X-ray pulsar XTE J1807-294, 



(N 

CN I focusing on its correlated behaviors in X-ray intensities, hardness ratios, pulse 

\^ [ profiles and power density spectra. The source was observed to have a serial of 

O I broad "puny" flares on a timescale of hours to days on the top of a decaying 

^ outburst in March 2003. In the flares, the spectra are softened and the pulse 

^ I proflles become more sinusoidal. The frequency of kilohertz quasi-periodic oscil- 

^ I lation (kHz QPO) is found to be positively related to the X-ray count rate in the 

^ I flares. These features observed in the flares could be due to the accreting flow 

inhomogeneities. It is noticed that the fractional pulse amplitude increases with 
^ the flare intensities in a range of ~ 2% — 14%, comparable to those observed in the 

' thermonuclear bursts of the millisecond X-ray pulsar XTE J1814-338, whereas 

it remains at about 6.5% in the normal state. Such a signiflcant variation of the 
pulse proflle in the "puny" flares may reflect the changes of physical parameters 
in the inner disk accretion region. Furthermore, we noticed an overall positive 
correlation between the kHz QPO frequency and the fractional pulse amplitude, 
which could be the flrst evidence representing that the neutron-star surface emis- 
sion properties are very sensitive to the disk flow inhomogeneities. This effect 
should be cautiously considered in the burst oscillation studies. 
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1. INTRODUCTION 

Neutron stars in low-mass X-ray binaries (LMXBs) accrete matter from a companion 
with a mass of less than 1 Mq via an accretion disk. In many models, the Keplerian accretion 
disk is supposed to be terminated at an inner radius Tin of a few Schwarzschild radii by 

e.g. relativistic effects, radiation drag, or neutron-star magnctosphere-disk interactions (see 
van der Klis 2004 for a review). The motion of the inner disk flow and the geometry of 
the innermost disk are still uncertain. The observations of the correlated behaviors in X- 
ray spectral and timing variabilities in LMXBs have provided the important probes of the 
accretion-flow dynamics and the stellar-disk interactions. 

Kilohertz quasi-periodic oscillations (kHz QPOs, in the frequency range of 200 Hz-1300 
Hz) have been observed in more than 20 accreting neutron star LMXBs (see van der Klis 
2000, 2004 for the recent reviews). Their frequencies are usually regarded to be associated 
with the Keplerian orbital frequency at some preferred radius related to rm. An important 
evidence for the movement of the inner edge of the disk under different mass accretion rate 
comes from the observation of a positive correlation of frequency vs. count rate on timescales 
of hours to days in some low- luminosity LMXBs, i.e. "atoll" sources (Mendez et al. 1999). A 
similar correlation (or in some cases an anti-correlation) has been also found in some high- 
luminosity LMXBs, i.e. "Z" sources, in the form of the correlations with curve length S^, 
along the track traced in an X-ray color-color diagram (Wijnands ct al. 1997; Yu, van der Klis 
& Jonker 2001; Homan ct al. 2002). The effect of stellar radiations on the disk movement 
has also been observed, e.g. in 4U 1608-52, where Yu & van der Klis (2002) found the 
anti-correlations between the kHz QPO frequency and the X-ray count rate associated with 
mHz QPO due to nuclear burning on the neutron star surface. 

Nearly coherent brightness oscillations have been discovered during thermonuclear X-ray 
bursts in some LMXBs (see Strohmayer & Bildsten 2003 for a review). Recent observations 
of the accreting miUisecond pulsars, SAX J1808.4-3658 (Chakrabarty et al. 2003) and XTE 
J1814-338 (Markwardt et al. 2003), confirm that the burst oscillation frequency is extremely 
close to the spin frequency. It supports the interpretation of the burst oscillation in terms 
of a "hot spot" on the stellar surface (Chakrabarty et al. 2003; Strohmayer & Bildsten 
2003), which triggers the efforts to investigate X-ray emission properties on the neutron 
star surface by folded pulse profiles (Strohmayer et al. 2003; Bhattacharyya et al. 2005; 
Watts, Strohmayer & Markwardt 2005). The best-fitting parameters of the pulsed fraction 
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and the first harmonic can not only determine the compactness of neutron star, but also 
help to figure out how the nuclear burning front (the hot spot) propagates on the stellar 
surface (Bhattacharyya et al. 2005; Watts, Strohmayer & Markwardt 2005). However, in the 
burst-oscillation studies, because of the rather limited photon statistics on the short burst 
timescale, many similar bursts have to be added to improve the statistics of the folded pulse 
profiles. The underlying problem is that these bursts may occur in different mass accretion 
rates. Till now, we know little about what effects of the accretion rate fluctuations can make 
on the neutron-star surface emissions or material depositions. 

This situation can be changed with the discovery of the accreting millisecond X-ray 
pulsars (MXPs, sec Wijnands 2005 for a latest review). Their detectable coherent pulsations 
allow us to get the pulse profile with good statistics over a long observational period when 
the energy spectrum (or hardness ratios) varies only slightly. According to the detailed 
spectroscopic and pulse-profile analysis on SAX J1808.4-3658, Poutanen & Gierhnski (2003) 
constructed a model for X-ray emissions from a hot spot on the surface/boundary layer 
of a rapidly rotating neutron star in the hard island state. They used it to constrain the 
compactness of the central star, the position of the emitting region and the size of the hot 
spot. Their investigations shed a light on further studies of the variations of the surface 
X-ray emission regions during the evolution of the accretion disk. 

Among the discovered 7 millisecond X-ray pulsars, XTE J1807-294 is the best candidate 
for us to investigate the impact of disk flow activities on the neutron-star surface emission 
for four reasons: (1) Besides SAX J1808. 4-3658, XTE J1807-294 is also a source which has 
been reported to have twin kHz QPOs ^ (Wijnands et al. 2003; Wijnands 2005; Linares 
et al. 2005). (2) The binary parameters of it have been calculated by Campana et al. 
(2003) and Kirsch et al. (2004) based on the XMM-Newton observation, so we are able to 
correct the photon arrival times for the orbital motions in producing pulse profiles. (3) 
The time-averaged energy spectrum of the source is found to be dominated by an optically- 
thin Comptonized component (Falanga et al. 2005), similar to that of SAX J1808. 4-3658 
(Gicrlinski, Done & Barret 2002), so that we can take the X-ray emission model of SAX 
J1808. 4-3658 as a direct reference. (4) There are no thermonuclear bursts reported in this 
source; thus we can focus on studying the effect of disk evolutions. In the left 6 MXPs, i.e., 
SAX J1808.4-3658, XTE J1751-305, XTE J0929-314, XTE J814-338, IGR J00291+5934 
and HETE J1900. 1-2455 (see Wijnands 2005 for a review), the above four conditions cannot 
be satisfied simultaneously. With the knowledge mentioned above, many basic features of 



^In completion of our paper, we noted that Linares et al. (2005) used the same data of XTE J1807-294 

as we analyzed. Our results are similar in kHz QPOs and different in the low frequency ranges, which might 
be due to the different choices of the data segments and different data grouping. 
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XTE J 1807-294 have been available as well, such as the shortest orbital period of ~ 40 
minutes and a relatively slow spin frequency of ~ 191 Hz (Markwardt et al. 2003). In 
addition, this source locates at 5°. 7 away from the Galactic center, with the best known 
position reported by Markwardt et al. (2003) based on a Chandra observation. Assumed 
a distance of 8 kpc, the source luminosity dropped from 1.3 x 10^''ergs~^ on February 28, 
2003, to 3.6 X 10^6 ergs"^ on March 22, 2003 (Falanga et al. 2005). 

We analyze the evolution of light curves, hardness ratios, pulse profiles and power den- 
sity spectra of XTE J1807-294, using Rossi X-ray Timing Explorer {RXTE) observations 
from February 27 to March 31, 2003. Our results firstly show that the positive kHz QPO 
frequency-count rate correlation on timcscalcs of hours to days is related to some "puny" 
flares, originated from disk flow inhomogeneities; the behaviors of the pulsed fraction and the 
flrst harmonic in the flares are different from those in the non-flares. However, in the whole 
investigated episode, there is a positive correlation between the pulsed fraction and the kHz 
QPO frequency, which could be the first evidence that neutron-star surface emissions are 
very sensitive to the disk flow inhomogeneities. Furthermore, we describe the RXTE obser- 
vations and data analysis method in Section 2; the results of timing analysis are presented 
in Section 3 and discussions of the results are shown in Section 4. Finally, we make a brief 
conclusion of our studies in Section 5. 



2. RXTE OBSERVATIONS AND DATA ANALYSIS 

We study the data obtained by the Proportional Count Array (PCA) on board RXTE 
during the March 2003 outburst decay of XTE J1807-294, using FTOOLS 5.3. 

To have an overview of the source evolution, we extract fight curves from the PCA 
standard 2 mode data (detector 2 only, all layers), and subtract the background produced 
by 'pcabackcst' version 3.0 with the "CM" faint-source model. We define soft color as the 
ratio of count rates in the 4.1 — 7.0 keV and 2.0 — 4.1 kcV bands and hard color as the 
ratio of count rates in the 13.3 — 30.2 keV and 7.0 — 13.3 keV bands. In selecting good time 
intervals for data analysis, we give a condition of PCU2_on.eq.l in addition to the standard 
conditions for faint sources. Then we extract the fight curves of PCU 0, 1, 3, 4 separately 
and compare them with the light curve of PCU2 respectively. Finally, we find the time 
intervals when PCU 4 meet instrumental fiares and delete such bad time intervals, making 
sure the following analysis are not infiuenced by instrumental problems. 

For Fourier analysis, we construct pulse profiles and power density spectra from the 
~ 0.95 /xs time-resolution PCA GoodXenon mode data (all available detectors, all layers). 
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correcting the photon arrival times for orbital motions of the pulsar and the spacecraft. 
When dealing with the GoodXenon mode data, we do not perform background subtraction 
or deadtime correction. 

For pulse-profile fitting we create folded light curves in 32 phase bins for every continuous 
event file. As a fit function of the pulse profiles we adopt a two-component Fourier function 
of the form 

= Co(l + ao sin[27r((/) - 0o)] + ai sin[27r(0 - 0i)/O.5]) (1) 

where Cq is the normalized count rate of the pulse profile; ao and Oi are the fractional 
amplitude of the fundamental and the first harmonic component, respectively; (J)o and 0i are 
the corresponding zero phase angles. We fit the pulse profiles and determine the errors of 
the fit parameters by using Ax^ = 2.7. We estimate the real averaged count rate from the 
2 — 60keV net light curve, produced by processing Standard 2 mode data (detector 2, all 
layers) within the corresponding intervals of every folded pulse profile. 

We construct power density spectrum (PDS) per observation using the data segments 
of 256 s and 1/4096 s time bins, so that the lowest available frequency is 1/256 Hz and the 
Nyquist frequency is 2048 Hz. We normalize the power spectra by the method of Miyamoto 
et al. (1991), and subtract the Possion white noise. In order to improve the statistics, we 
add the observations into groups when these observations are adjacent in time and the power 
spectra remain the same. The division of the groups can be found in Table 1 and are depicted 
in the top panel of Fig. 1. Most of our power spectra show a power-law-like component in 
mHz-Hz frequency range, which could be the very low frequency noise (VLFN). Furthermore, 
in our PDS, some additional weak noises can be measured with good statistics as shown in 
Section 3.3. 

We use a power law plus multi-Lorentzian function (Belloni, Psaltis & van der Klis 
2002; van Straaten et al. 2002) to fit the spectra. Before fitting the power spectra, we remove 
frequency bins containing the pulse spike. To visualize the characteristic frequencies (fmax) of 
noises or QPOs, we plot the power spectra in the power times frequency representation {lyPu), 
where P,^ is the power spectral density and u is the Fourier frequency. In this representation, 
the Lorentzian's maximum occurs at z/max (^'max = \/ J^o + where z/max is the characteristic 
frequency, z/q is the centroid frequency and A is the half width and half magnitude (HWHM) 
of the Lorentzian). We represent the Lorentzian relative width by Q defined as 1/0/2 A, and 
the root- mean-square fractional amplitude by rms. 
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3. RESULTS 

3.1. Broad "Puny" Flares in 2003 March 

It has been found that the light curve of XTE J1807-294 drops exponentially from 
February to the end of March of 2003 (Falanga et al. 2005) and then becomes flattened 
instead of being much steepened as in SAX J1808.4-3658 (Gilfanov et al. 1998; Gierlinski, 
Done & Barret 2002) and XTE J1751-305 (Gierlinski & Poutancn 2005). Figiu^c 1 shows 
that in the former exponentially decaying episode, at least four intensity fluctuations occur, 
each lasting beyond several thousands of seconds, with small X-ray count rate variations. 
To make quantitative analysis, we flt the light curve by an exponential plus multi-Gaussian 
model. The full width and half magnitude (FWHM) of each Gaussian and the ratio of 
the magnitude of the Gaussian to the exponential value at the same time are (0.53 day, 
7.18/45.42), (0.95 day, 16.8/35.22), (1.67 day, 7.53/20.75), (0.71 day, 12.33/17.16), respec- 
tively. The features of long duration and low amplitude make these fluctuations different 
from type I thermonuclear bursts and super nuclear outburst (Strohmayer & Brown 2002; 
Strohmayer & Markwardt 2002) on the neutron star surface. We thus call them "puny" 
flares for convenience. In the following, it can be noticed that these "puny" flares are also 
featured by the enhanced soft X-ray emissions and increased pulsed fractions. 

The second panel of Fig. 1 shows that the soft color drops in the flares, especially at 
the top of the largest flare of Mach 14-15. By comparing the light curves in the four energy 
bands (2.0 - 4.1 keV, 4.1 - 7.0keV, 7.0 - 13.3 keV and 13.3 - 30.2 keV), we flnd that this 
is due to the stronger soft X-ray emissions in the band of 2.0 — 4.1 keV. The color- intensity 
diagram of Fig. 3 indicates that, except of the "puny" flares, the soft color decreases with 
the decay of the outburst, whereas the hard color changes a little in the whole investigated 
episode. 

We can deduce from the almost invariant hard color that the hard spectral component 
changes only slightly in the whole episode. The fltting results of the combined spectra of 
XMM-Newton, RXTE and INTEGRAL observations (Falanga et al. 2005) present that hard 
emissions from the hot and optically thin Gomptonized plasma dominate the energy spectra; 
the contribution of thermal emissions from the accretion disk at February 28 and March 22 
are so small that they are almost negligible in PGA analysis (above 2.5 keV). Therefore, 
the PGA X-ray emissions can be considered to be mainly generated from the neutron star 
surface/boundary layer. 
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3.2. Variations of Pulse Profiles 

We fit 86 pulse profiles with a two-component sine function described in Section 2. 
The fitting is good in most cases, with x^/dof ranged in 32/59 — 80/59. The situations of 
X^/dof > 1.3 only occur in three pulse profiles. We plot the evolutions of the best-fitting 
parameters of Qq and ai in the fourth and the fifth panel of Fig. 1, respectively. It can be 
seen that oq is a rather unambiguous indicator of the "puny" fiares than the X-ray intensities 
and soft colors, not only because of its small errors, but also because oq increases significantly 
at the top of the fiares despite how weak the fiare is. 

Figure 2 shows the different behaviors of Oq and Oi with respect to the averaged count 
rate of the pulse profile. Combining the fourth panel of Fig. 1, we find that in the normal 
state, the first harmonic is relatively strong, but it becomes weaker (with ai decreasing from 
1.4% to about 0.9%) in the decay of the outburst, whereas oq keeps at ~ 6.5%; in the fiares, 
ao is positively correlated with the averaged count rate, whereas Oi is generally in a lower 
level around 0.7% over the fiares. The variational range of Oq is ~ 2% — 14%. 

3.3. Aperiodic Timing Behavior 

We fit the power density spectra by one power-law component plus three to seven 
Lorentzian components. Errors on the fit parameters are determined using Ax^ = 1. We 
list the characteristic parameters of the Lorentzians, e.g. Umax, Q and rms in Table 2, the 
power-law index a of the VLFN and the x^/dof of the fitting in Table 3. In our investigated 
episode, kHz QPOs have been detected in groups 1-7 and 9-10, and undetectable in the other 
groups on account of the poor statistics of the power spectra above ~10 Hz. The typical 
power spectra of XTE J1807-294 are shown in Fig. 4, where group 1 is chosen to stand for 
the non-fiare state and group 5 is taken as a representation of the fiare state. 

We plot the characteristic frequencies of the Lorentzians versus the upper kHz QPO 
frequency in Fig. 5. It shows that XTE J1807-294 has a frequency-frequency correlation 
similar to that of SAX J1808.4-3658 (van Straaten, van der Khs & Wijnands 2005). Most of 
the Lorentzian components can then be identified in the scheme of Belloni, Psaltis & van der 
Klis (2002) and van Straaten, van der Klis & Mendez (2003). In Fig. 5, L„ is the upper kHz 
QPO; L^ is the lower kHz QPO; LhHz is the hectohertz Lorentzian, a broad Lorentzian with 
a frequency around 150 Hz; L;, is the break Lorentzian, a band-limited noise component; Lh 
is the low frequency Lorentzian just above L5. It should be noted that when the frequency 
of Li approaches the range of 100 ~ 150 Hz, it becomes difficult to be distinguished from 
^hHz, especially when one of them is not detected in the spectra, e.g. groups 3 and 4 (see 
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Fig. 5). In Table 2 we put '?' on their identifications. 

In Fig. 5 we multiply the frequencies of twin kHz QPOs of Z sources GX 5—1 and GX 
17+2 by a factor of 1/1.5 for comparison. Their shifted frequencies are almost consistent 
with the frequency-frequency distribution of the twin kHz QPOs of XTE J1807-294. We 
draw Q versus fmax for both and in Fig. 6. It proves that the frequency shift exists 
in the upper kHz QPO frequencies in XTE J1807-294 as in the case of SAX J1808.4-3658 
(van Straaten, van der Klis & Wijnands 2005). 

Except of the five major Lorentizans of L^, L^, LhHz, and L^, the power spectra of 
XTE J1807-294 present some additional components (Fig. 5). One of the components is 
the Lorentizan of Lh2, with frequency nearly double of Lh- If Lh corresponds to the HBO 
component of Z sources (see van der Klis 2004 for a review), L/j2 would correspond to the 
HBO harmonic. Such harmonic has never been reported in atoll sources or other accreting 
millisecond pulsars. Another additional component is the Lorentizan of Llfn, which takes 
the position of the shifted low- frequency noise (LFN) component of GX17 + 2 in the diagram 
(Fig. 5). The fitted power law index of ulfn vs. is 3.42 ±0.29, larger than the index of Uf, 
vs. Uu of ~2.7. We use 'Ftest' to check the probabihty of the components whose errors are 
less than 2a. All of the weak noise components hsted in Table. 2 have Ftest probabilities 
less than 0.03. 

The simultaneous detection of and L^^ in XTE J1807-294 makes it evident that 
Llfat is a different component from Lj,. We infer from this finding that the detected LFN 
in GX 17 + 2 (Homan et al. 2002) should be a different component from those observed in 
other Z sources. The shift factor of about 1.5 in Llfn as well as in L„ and L^ imply that 
these three components could originate from the same inner disk region. 

3.4. Correlations Between Count Rate, Pulse Profile and kHz QPO 

We draw in Fig. 7 the upper kHz QPO frequency (u^) versus the averaged count rate 
and the fractional pulse amplitude (ao) of the group, respectively. The averaged count rate 
and oo of every group were calculated using the same methods described in section 3.2. 

We can see from Fig. 7 and Fig. 5 that, while the frequencies of the Lorentizans decrease 
with the count rate in the non-flare state of groups 1-3 and increase with flare intensities in 
groups 4-6 and 7-10, they do not vary monotonoTisly with the X-ray count rate in the whole 
episode. The positive frequency-count rate correlations in the rise of the flares seem to be 
steeper than that measured in the normal outburst decaying state. However, in the whole 
investigated episode there exists a noticeable positive u-^ vs. ao relation with a linear slop 
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coefficient of ~ 35 (Fig. 7). 

4. DISCUSSION 

Our motivation to do timing analysis on XTE J1807-294 is to probe the effects on 
tfie neutron-star surface X-ray emissions made by the disk flow activities, so we focus on 
analyzing the variations of pulse profiles and power density spectra during the evolution. 

We find several pairs of twin kHz QPOs in the power density spectra of XTE J1807-294. 
The frequency separation Ai/ ranges from ~ 179 Hz to ~ 247 Hz, close to the spin frequency 
of Vs ~ 191 Hz. Several kHz QPO models, considering the coupling of the neutron-star spin 
to the Keplerian orbit motion of material in the accretion disk, have been put forward to 
explain this specific feature. We notice that most of these models concern inhomogeneous 
accretion fiow by for example accreting blobs, density fluctuations in the innermost part 
of the accretion disk, or high-density loops confined along a magnetic field line. In the 
modified sonic-point beat-frequency model (Lamb & Miller 2001, 2003), the accreting blobs 
are introduced to explain the deviation of Az/ from the constant value of the spin frequency 
by considering their relative motions to the neutron star. In the resonant disk oscillation 
models (Kluzniak et al. 2004; Lee, Abramowicz & Kluzniak 2004), the twin kHz QPOs 
are explained by a nonlinear resonance excited by coupling disk oscillation modes to the 
neutron star spin. Petri (2005) proved further that a rotating misaligned magnetic field of 
a neutron star or a rotating non-axially symmetric gravitational potential can provide the 
coupling; their perturbations on the accretion disk can give rise to very pronounced density 
fiuctuations at the inner edge of the disk. Most recently, Li & Zhang (2005) pointed out 
an alternative interpretation for the lower kHz QPO in terms of the standing kink modes of 
magnetic loop oscillations at the inner edge of the disk, where the loops with high-density 
plasma and small cross section can be produced from the reconnection of the azimuthal 
magnetic field lines. 

The variation of kHz QPO frequencies can be explained by the movement of the inner 
edge of the accretion disk in the radiative disk truncation model (Miller, Lamb & Psaltis 1998; 
van der Klis 2001). For example, the observed positive frequency-count rate correlations on 
timescales of hours to days in some atoll sources have been interpreted as the changes in the 
inner disk radius due to disk accretion rate fiuctuations. We find a similar frequency-count 
rate correlation in XTE J1807-294, but we notice that the shift of the QPO frequency is 
related to some broad "puny" flares. The observed frequency- count rate correlation in the 
rise of the flares seems to be steeper than in the normal state. This result means that (1) 
in the case of the same count rate, the kHz QPO frequency measured in the rise of the 
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flare is larger than that in the normal state; (2) for a certain kHz QPO frequency range, 
the variation of the count rate is smaller in the flare than in the normal state. A possible 
explanation is that there are blobs of high mass densities at the inner edge of the accretion 
disk and the they can persist in time scales of hours to days. The formation and the detailed 
structure of such a kind of blob are still uncertain. We refer to the hydrodynamic simulation 
result by Petri (2005) mentioned in the above paragraph for instance. It might be expected 
that the relative intensity of the flares to the nearby normal state could be a measure of the 
mass density of the blob; while the limited time intervals of the flares could be related to 
the size of the blobs. However, besides the variation of mass accretion rate due to accreting 
the blobs, there are other mechanisms which can influence the X-ray intensity and kHz QPO 
frequency, such as nuclear burning on the neutron star surface (Bildsten 1993; Yu & van der 
Klis 2002). To make a further quantitative investigation of the blobs from the observation, 
we have to distinguish their different effects firstly. 

Accreting inhomogeneous flow as one of the possible origins of flares has been suggested 
by Moon, Eikenberry & Wasserman (2003) for LMC X-4. In the "broad flare" of LMC X-4, 
the spectrum is softened and the pulse proflle becomes simple sinusoids. These features are 
also observed in the broad "puny" flare of XTE J1807-294. In the massive X-ray binary LMC 
X-4, the X-ray emission mechanisms of the magnetic neutron star (with surface magnetic 
field strength B ~ 10^^ G) are too complex to be modelled. However, the spectrum of XTE 
J1807-294 can be well fitted by an absorbed black-body plus Comptonization model; thus 
we can hope to learn more about the X-ray emission properties by fitting the pulse profiles of 
XTE J1807-294 as done by Gierhnski, Done & Barret (2002); Gierhnski & Poutanen (2005) 
for SAX J1808.4-3658 and XTE J1751-305. 

Our preliminary results of the pulse-profile analysis present that both oq and Oi be- 
have differently in the fiares and in the non-fiares. As for the latter case, while the X-ray 
count rate decreases, Oq keeps at about 6.5%, whereas Oi decreases from ~ 1.4% to ~ 0.9%. 
However, in the former case, ai is generally in a low level of 0.7%, whereas ao increases 
with the flare intensities. The value of ag spans in the range of 2% — 14%, compara- 
ble to those observed in some thermonuclear bursts of XTE J1814-338 (Strohmayer et al. 
2003; Watts, Strohmayer & Markwardt 2005) and some non-pulsing LMXBs (Muno, Ozel & 
Chakrabarty 2002). According to the X-ray emission model of Poutanen & Gierhnski (2003), 
the variation of the pulse profile can be determined by the changes of the surface-emission 
parameters, such as the hot spot size, position, and even emission patterns, e.g. the optical 
depths of the black-body emission and the Comptonization, respectively. We will put the 
detailed spectroscopic and pulse-profile modelling analysis in a future work. According to 
this model, the above results have already shown that the accretion emissions at the neutron 
star surf ace/boundary- layer are very sensitive to the accretion fiow inhomogeneities. The 
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most direct evidence is the observed positive correlation between kHz QPO frequency and 
Qq in the investigated episode. 

Although we infer from the above results that the broad "puny" flares could be mainly 
due to accreting blobs, we cannot exclude the possibility of the existence of some special 
mode of nuclear burning on the neutron star surface. Bildsten (1993) firstly proposed that a 
time-dependent nuclear burning (slow fires) on patches of the neutron star surface can occur 
in an intermediate accretion regime, where the accretion rate is sub-Eddington but is still 
too high to allow type I thermonuclear burst, e.g. 5 x 10^^° MqJT^-^ < M < 10~® MqJt'-^ 
for a pure helium burning case (Bildsten 1995). For XTE J1807-294, the conditions for slow 
fires seem to be satisfied qualitatively: no type 1 thermonuclear burst has been observed and 
the source luminosity drops from 1.3 x lO^^ergs"^ to about 3.6 x lO^^ergs"^ (assuming the 
distance of 8 kpc) in the investigated episode (Falanga et al. 2005). However, the property 
of the companion star is still uncertain. It could be a He dwarf with a low inclination angle 
(i < 30°) located at ~ SkPc, or a C/0 dwarf with high inclination angle (60° < i < 83°) at 
about 3 kPc (Falanga et al. 2005), which makes it difficult to give a quantitative estimation of 
the ignition and the spreading of the fires on XTE J1807-294. Even though, the observations 
of the low-level luminosity variations on timescales of 10'^ s and the existence of the very low 
frequency noise (VLFN) in the mHz-Hz range are all consistent with the features of the slow 
fires predicted by the theoretical investigations of Bildsten (1995). If the slow fires exist, the 
position and the size of the black-body radiation on the neutron star surface will certainly be 
influenced, making impacts on the soft-band pulse profiles and time lags. We will consider 
this effect in another work. 



5. SUMMARY AND CONCLUSION 

XTE J 1807-294 shows significant variations of pulse profiles and kHz QPO frequencies 
during the source evolution in March, 2003. It thus provides a good chance to probe the 
neutron-star surface/boundary layer X-ray emissions, the inner disk fiow activities and their 
correlations. 

The main results we get from RXTE data analysis are related to some specific "puny" 
fiares, featured by low-amplitude intensity fiuctuations, hours-to-days long timescale, stronger 
soft emissions and significant increase of fractional pulse amplitude with the fiare intensities. 
A positive kHz QPO frequency-count rate correlation is also found to be connected with the 
flares: the correlation seems to be steeper in the rise of the flares than in the normal state. 
However, in the whole investigated episode, the kHz QPO frequency is positively correlated 
with the fractional pulse amplitude. In the normal state, the fractional pulse amplitude 
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keeps constant and the first harmonic content decreases with the outburst decay. 

We propose that accreting high-mass-density blobs could be responsible for all of the 
above features of the "puny" flares, and conclude that neutron-star surface emissions should 
be very sensitive to the disk flow inhomogeneities in the hard state. Since the variation of 
the fractional pulse amplitude in the flare is comparable to that in the thermonuclear burst, 
the effect of the accretion flow inhomogcncity on material depositions at the neutron star 
surface should be carefully considered in the study of thermonuclear burst oscillations. 

Detailed spectroscopic and pulse- profile analysis on XTE J 1807-294 will be done to 
determine the physical parameters of the surface/boundary-layer emissions in the whole 
investigated episode in a subsequent work. 

We would like to thank the anonymous referee for his or her very helpful comments, 
and thanks are also due to M. Linares and S. Zhang for their helpful discussions on the 
RXTE data analysis. This research has made use of the data obtained from the High 
Energy Astrophysics Science Archive Research Center provided by NASA's Goddard Space 
Flight Center. This work is subsidized by the Special Funds for Major State Basic Research 
Projects and by the National Natural Science Foundation of China. 
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Fig. 1. — Evolution of the light curve, soft/hard color, fractional amplitude of the pulsation 
(oo) and the first harmonic component (ai) of XTE J1807-294 in March, 2003. The hght 
curve of 2. — 30.2 keV is in unit of counts s~^ PCU~^ (with time bins of 256 s); the soft color is 
defined as i"[4.i-7.okeV]/%.o-4.ikeV], and the hard color is defined as i"[i3.3-30.2keV]/-^[7.o-i3.3keV] 
(with time bins of 16482 s). The horizontal coordinate stands for the date of the observations, 
where -1 and denote the days of February 27 and 28, respectively. The four arrows in the 
fourth panel indicate the positions of the flares. 




Fig. 2. — Best-fitting parameters of the fractional pulse amplitude Gq (up panel) and the 
fractional first harmonic amplitude oi (bottom panel) with respect to the averaged net count 
rate of the pulse profiles. 
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Fig. 3. — Hardness intensity diagram of XTE J 1807-294. The group numbers for every 
power density spectra are noted in the top paneL 
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Fig. 4. — Power density spectra in vPy of group 1 (left) and group 5 (right) for XTE J1807- 
294. Solid lines mark the best-fit model, and dotted lines are the components of every 
Lorentzian. The components are identified by their distributions in the frequency-frequency 
diagram of Fig. 5. 
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Fig. 5. — Characteristic frequencies of the power spectral components of XTE J1807-294 
versus the upper kHz QPO frequency. The data points of SAX J1808. 4-3658 (van Straaten, 
van der Khs & Wijnands 2005) are drawn for comparisons. Also shown are of the Z sources 
GX 5-1 (Jonker et al. 2002) and GX 17 + 2 (Homan et al. 2002), and I^lfn of GX17 + 2, 
where the frequencies of these Z sources are all shifted by a factor of 1/1.5. 
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Fig. 6. — Q value of L„ (left) and Lh (right) of XTE J1807-294 (triangles) plotted versus 
the corresponding characteristic frequencies respectively. Data points of SAX J1808. 4-3658 
(van Straaten, van der Khs & Wijnands 2005) are plotted by grey squares for comparisons. 
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Fig. 7. — Characteristic frequencies of the upper kHz QPOs versus the averaged count rate 
(left) and the fractional pulse amphtude Qq (right) of groups 1-7, 9-10. The group numbers 
are noted in the right panel. 
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Table 1: Observational parameters of the 16 groups in Fig. 1. 



No. Observation ID 


Start 


End 


Exposure (s) 


1 70134-09-02-00, 


Feb. 27 11:44:32 


March 1 00:44:08 


9.712E+03 


70134-09-02-01, 








80145-01-01-01 









(80145-01-) 



2 


01- 


-03, 


01- 


■04, 01-00 


March 1 


12:44:04 


March 2 


03:43:44 


1 


133E+04 


3 


02- 


-00 






March 4 


22:53:04 


March 5 


01:31:44 


6 


032E+03 


4 


02- 


-03, 


02- 


-02, 02-01, 


March 6 


16:04:00 


March 8 


00:30:40 


1 


912E+04 




03- 


-03, 


03- 


-02 














5 


03- 


-00 






March 8 


23:42:08 


March 9 


03:19:44 


8 


096E+03 


6 


02- 


-04, 


02- 


-06, 01-05, 


MarchlO 


00:28:32 


Marchl3 


01:35:44 


4 


984E+04 




03- 


-01 


















7 


02- 


■05, 


04- 


■01 


Marchl3 


18:27:12 


Marchl4 


01:44:32 


1 


536E+04 


8 


04- 


■00 






Marchl4 


21:17:20 


Marchl5 


01:55:44 


9 


696E+03 


9 


04- 


-02, 


04- 


■10 


Marchl6 


00:06:08 


Marchl6 


04:15:44 


7 


792E+03 


10 


04- 


-03 






March 16 


20:37:20 


Marchl7 


00:46:40 


9 


696E+03 


11 


04- 


-04, 


04- 


-09, 04-11, 


Marchl7 


18:45:20 


Marchl9 


22:14:40 


1 


598E+04 




04- 


-05, 


04- 


-07 














12 


04- 


-08, 


05- 


-01, 05-03, 


March20 


20:51:28 


March23 


19:26:40 


2 


941E+04 




05- 


■02, 


05- 


-00, 05-11, 
















05- 


■05 


















13 


05- 


■06, 


05- 


■04, 05-08 


March25 


00:46:24 


March27 


00:43:44 


6 


128E+03 


14 


05- 


-07, 


05- 


-10, 05-09, 


March27 


16:03:28 


March29 


19:29:36 


1 


984E+04 




06- 


-00, 


06- 


-07, 06-01, 
















06- 


-10 


















15 


06- 


-02 






March30 


14:35:28 


March30 


16:56:32 


5 


840E+03 


16 


06- 


-03 






March31 


11:48:32 


March31 


17:20:32 


9 


664E-f03 
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Table 2. Characteristic parameters of the fitted Lorentzians of groups 1-10. 



No. 






Li:,FN(L62) 


Lb 


Lfcl 




Lfc2 








1 


"^mBx{Hz) 
Q 

rms(%) 


- 


, q, + 0.45 

0.72l°;« 
4.5ll;3 


6 60+°'''l 
°-°"-0.67 

„ 0^+0.23 

8.0 ± 1.0 


- 


24.7 ± 1.0 
1.37 ± 0.25 
11.2+3-3 


64+1* 

"■76l?,:S 

10.6l^2-« 


- 


240^33 

, „_+0.67 
^■'"-0.45 
10.6i3-° 


460.4 ± 2.5 
6.9 ± 0.7 
13.5±1.1 


2 


Q 

rms{%) 


- 




4.66 ± 0.61 
0.42 ± 0.11 
8.9 ± 1.1 


- 


17 41+0.80 

l-°'-0.35 
8.1 ± 1.7 


42 2+"^^ 
^''■^-lO.S 

0^44t°o:?o 
14.49l3;« 


- 


200+ 

2-i9iJ:fo 

8-0913-];* 


378.8 ± 2.4 
,S.74 ± 0.56 
12.8 ± 1.0 


3 


i'max(Hz) 

Q 

rma(%) 


- 


- 


3 20+"-''^ 
0.64l°;3? 
6.4ti-*3 


- 


17.8l-« 

0.47l°;ll 

14.0 ± 1.9 




122+*" 

o.98±S:;!^ 

9.2i-^-^ 


> ? 

> ? 

> ? 


322 ± 11 

q qQ+1.31 
O. jy_Q gg 


4 


i'max(Hz) 

Q 

rma(%) 


- 


- 


5 64+" '*3 

49+° l'5 
"•'"'_0.12 

1.9 ± 1.3 


14.43lli:^,^ 


18.0+;-^ 

2 gi + 1.60 
^-Oi_0.79 
B.6 ± 1.7 


30.8 ± 7.6 
0.68ig:fg 

io.9j:|;3 


162+*^ 

70+°^*8 
"■ '"-0.29 

ii.sjil;! 


> ? 

> ? 

> ? 


368.9 ± 6.6 
, q<,+0.81 
•3-""-0.64 
11.8 ± 1.5 


6 


i^max(Hz) 

Q 

rm,s(%) 


1.68 ± 0.06 
4 7+3-1 

"'■"-0.7 


3-°O-0.23 

3 6+3 '' 
•^■°-1.8 


11 9+^''' 
11 4+1-8 




39.2+;J-3 

2 24+0-133 
''-■^*-0.46 

"-"-1.4 




0.69l°;f, 

9.8iti 


360 ± 11 

4 6+1^* 
^■"-1.2 

9 1+23 
"■^-2.1 


545.3 ± 2.1 
13.5+2'3 
10.4+1^? 


6 


l^max{Hz) 

Q 

rms(%) 






6-87+093 
r. ,0 + 0.11 
" -5»_o.09 
9.9 ± 1.1 




27.3 ± 2.0 
90+"-^^ 
11.6 ± 1.3 






237+°° 

57+"^31 
"■='-0.26 

13.5_2'2 


458.2+*;^ 
6.6 ± 1.1 
10 1 + 1^* 


7 


i^max(Hz) 

Q 

rms(%) 






11 6 + ^-^ 
0.38 ± 0.12 

V2.ii+_\l 






41-5i'2-3 
, __+0.i30 
'-■ ' ' -0.43 
9.2+l,^? 




S4n+82 

■J^"-65 

i.26ij:°i 

10.8it.°3 


587.0+^-^ 
9.413;? 


8 


i'max(Hz) 

Q 

rms(%) 






10.8+J-^ 

"-■="-0.07 
14 1 + 1-^ 






^■4it.l 








9 


i'max(Hz) 

Q 

rms(%) 






9.9+^-^ 

40+'''^^ 
"-^^_0.19 

10 7+=* 
^"-'-1.9 




39-0+^2 
"-"-2.3 






275+»^ 

"■ ' " — 0.36 
14 Q+^-^ 


501.4+«^5 

io.8±«3:«g 
8.siti 


10 


i'max(Hz) 

Q 

rms(%) 






4 97+»-«« 
—0.82 

0.00 (fixed) 
12.9^1;? 




19.2+^-^ 

99+°-*= 
"■"" — 0.30 

9 2+2^1 
9^'=-1.9 








368+52 

98+°^*2 
"-"° — 0.31 

16.6111 



Note. — In groups 3—4, ' > ?' is used in the case where can hardly be distinguished from I^^Hz- 



Table 3. Fitted power law index of the VLFN component together with and degree of 



the freedom (d.o.f) of the fitting. 



No. 


1 


2 


3 


4 


6 


OL 


-2.04 ± 0.16 


-1.99 ± 0.08 


-1.78 ± 0.09 


-1.86 ± 0.14 


-1.74 ± 0.18 




216 


214 


208 


408 


221 


d.o.f 


189 


204 


195 


232 


229 


No. 


6 


7 


8 


9 


10 


a 


-1.90 ± 0.06 


-1.73 ± 0.09 


-1.95 ± 0.21 


-2 34+'J^3li 
^■■^4-0.41 


-1.94 ± 0.15 




397 


289 


261 


343 


267 


d.o.f 


228 


238 


226 


203 


239 



